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We present a detailed study to show the possibihty of approaching the quantum ground-state 
- - - of a hybrid optomechanical quantum device formed by a Bose-Einstein condensate (BEC) confined 

CO , inside a high-finesse optical cavity with an oscillatory end mirror. Cooling is achieved using two 

. experimentally realizable schemes: back-action cooling and cold damping quantum feedback cooling. 

»' ' In both the schemes, we found that increasing the two body atom-atom interaction brings the 

mechanical oscillator to its quantum ground state. It has been observed that back-action cooling is 
C\ . more effective in the good cavity limit while the cold damping cooling scheme is more relevant in 

D ' the bad cavity limit. It is also shown that in the cold damping scheme, the device is more efficient 

[JL^ ' in the presence of BEC than in the absence of BEC. 
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I. INTRODUCTION 



' Optomechanical cooling of micro- and nano-mechanical resonators to their quantum ground state, is used in a wide 
, variety of sensitive measurements such as detection of weak forces [l], |2| , small masses Q and small displacements (3] ■ 
^ ; Such mechanical resonators can also be used as tools for quantum metrology 0, Q or as a medium to couple hybrid 
quantum systems 0, II]- In recent years, different types of optomechanical systems have established considerable 
> cooling of the vibrational modes of mechanical resonators interacting with optical cavity [1, In such systems, the 

a[ cavity exerts radiation pressure on the mirror. The dynamics of the mirror is coherently controlled by an external 
pump laser. Therefore, by modifying the pump power, a strong coupling regime can be easily acheived In such 
' a strong coupling regime, the mechanical oscillator can be cooled to its quantum ground state using the dynamics 
(~| , of back-action (12i . ilS] . Experimentally, significant cooling of the mechanical mode of the resonator coupled to an 
O ' optical cavity has been attained by using two different ways of radiation pressure interaction between the intracavity 
, ' field and a vibrational mode : back-action cooling [3 Ell and cold damping feedback [l^, HI, [13 ■ Self cooling 
of the mechanical resonator via dynamical back-action arises due to its interaction with the optical cavity through 
7-H ' radiation pressure @, HO- Depending upon the laser detuning, the correlations, induced by the cavity delay, 
^ . between the Brownian motion of the resonator and the radiation pressure force lead to either cooling or amplification. 
On ■ Experimentally, a single mechanical mode has been cooled using these effects (T8l - [20j . With the back-action cooling 
. or self cooling, there is self modification in the dynamics of the mechanical s yste m as the off resonant operation of 
' the cavity results in the retarded dynamical back-action on the system IT8| - |2T| . Cold damping quantum feedback 
technique is used to cool the oscillating mirror by applying a viscous force on the mirror via radiation pressure generated 



by another intensity modulated laser beam on the back of the mirror [23 - l25| . This technique increases the damping of 
the system without any additional noise |2Q . i27\ . Other feedback scheme requires very large mechanical quality factor 
[28j for the ground state cooling of the oscillator whereas cooling can be much more conveniently achieved in the 
cold damping scheme. The optomechanical systems composed of ultracold atoms enclose within a cavity has grabbed 
much attention[29~Sl|- The resonance frequency of the cavity is altered due to the strong interaction of the condensed 



k> atoms with the cavity mode [42|. An all-optical transistor based on a coupled Bose-Einstein condensate (BEC) cavity 
}_] . system [3l|, has been proposed [l^l ■ An optomechanical system consisting of BEC in an optical cavity has been 
■ investigated recently By virtue of interaction of BEC with mechanical oscillator, such hybrid system was helpful 
in acheiving state engineering of the mechanical mode of oscillator, which can be easily controlled and highy insensitive 
to noise effects. Experimentally, interaction between the ultracold atoms and vibrating membrane has been studied 
(46j . In our previous work, we have investigated how the stochastic cooling feedback technique together with a gas of 
condensate confined in an optical cavity can be used to detect weak forces and coherently controls the sensitivity of 
this hybrid optomechanical system [47| . 

Motivated by these interesting features in the field of cavity opto-mechanics and BEC, we compare the two schemes 
i.e., back-action cooling and cold damping quantum feedback, for cooling an optomechanical system in which a BEC 
is coupled to an optical cavity with a movable mirror. By comparing these two schemes for such a system, we found 
that the cold damping quantum feedback scheme cools the mechanical oscillator to its ground state in the bad cavity 
limit {k » ujm) i-e, when cavity bandwidth is greater than the mechanical frequency while the back-action cooling 
is efficient in the opposite limit of good cavity (/t << a;„i) i.e., when cavity bandwidth is smaller than the mechanical 
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frequency. We also show that for both the schemes, in the presence of atom - atom coupling, the interaction of the 
vibrational mode of the mirror, single mode of the intracavity field and the condensate fluctuations give rise to normal 
mode splitting. 



II. MODEL HAMILTONIAN 



In this section, we introduce our model and describe the Hamiltonian for our system. Our model consists of an 
elongated cigar shaped gas of N two- level BEG atoms of ^'^Rb in the \F — 1 > state with mass m and transition 
frequency Wa of the \F = 1 > \F' = 2 > transition of the D2 line. The BEG atoms strongly interact with a 
single one-dimensional quantized cavity mode of frequency ujc- The cavity has one mirror fixed and another mirror 
movable of mass Mm which acts as a mechanical oscillator that is free to oscillate at mechanical frequency Wm- The 
coherent laser field with frequency Up drives the system through the fixed cavity mirror with amplitude 77'. The 
cavity and the laser field are treated quantum mechanically as the standing wave and the pump laser are weak (small 
photon number) . The single-particle Hamiltonian in the rotating wave and dipole approximation for the system under 
consideration is read as ,291: 



Ho — Ha + He + Hm + Hca + Hem, (1) 

where Ha — ^ represents the kinetic energy of the condensate. He is given by —HAcb^b — ihrj (6 — b^) in which 
the first term gives the energy of the light mode, with lowering (raising) operator b (6^) ([6, 6^] — 1) and cavity-pump 
detuning Ac{— LOp — lOc)- The second term in He is the energy due to external pump laser. A band pass filter is 
used in the detection scheme [48!] such that only a single vibrational mode can be considered by neglecting the several 
other mechanical degrees of freedom arising from the radiation pressure. Hfyi — ^CUrn c* c gives the energy of the single 
vibrational mode of the cantilever with annihilation (creation) operator c (c^) ([c, c^] = 1). As a result of the pressure 
exerted by the intracavity photons on the mirrors, there is an optomechanical coupling between the vibrating mirror 
and the cavity field. Hence a force is exerted on the movable mirror by the electromagnetic field which depends on 
the number of photons in the cavity. This field is phase shifted by 2klm where k is the wave vector of light field and 
Im is the displacement of the mirror from its equilibrium position. The scattering of the photons into other modes 
of the electromagnetic field can be ignored and hence validating the single mode approximation of the cavity field. 
Also, a reservoir having finite temperature T is attached to the movable mirror. In the absence of the radiation 
pressure exerted by the optical cavity field, the cantilever can undergo Brownian motion due to its contact with 
thermal environment. Hca = cos^(A: x){hvQb% -f Vd) describes the interaction of the light field with the BEG which 
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involves the classical potential Vd, where vq is the optical lattice barrier height per photon given by uq = with 
go and Aa(= LOp — tUa) representing the mirror- atom coupling and atom-pump detuning respectively. We take f o > 
as in this case condensate atom moves towards the nodes of the optical field due to which the lowest bound state is 
localized at these positions. Therefore compared to the case uq < 0, the interaction between the optical cavity field 
and the condensate atom is decreased for -uq > 0. Hem — —heuJmb^b{c + c^) represents the energy due to non linear 
dispersive coupling between the position quadrature of the movable mirror and the intensity of light field where e is 
the mirror-photon coupling. The Hamiltonian in second quantized form with the two-body interaction is given as : 



H = J dr^\r)Ho^{r) + ^^^^^ J df ¥ {r)-^'< {r^^ (r)'^ {f) , (2) 

where ^'(r) is the atom field operator and as is the s-wave scattering length. Using 'i'{r) — aiw{r — fl), we derive 
the Bose-Hubbard Hamiltonian. Here, w{r — fl) is the wannier function and aj is the corresponding annihilation 
operator for the ith bosonic atom. Experimentally, the optical lattice depth is tuned such that time scales over which 
the experiment is performed is smaller than the time scales over which tunneling takes place, hence the tunneling of 
the atoms into neighbouring wells can be neglected in deriving the Hamiltonian. Therefore, Hamiltonian reads as : 

H = Kg ajQj — hAJj'b + fkdmc'c + Po{hvob^b + Vd) Q^aj — heuJmb''b{c + c^) 

j 3 

- ih-q [b-b^) +^-^'^a]a]aja.j, (3) 

'j 
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where Kq = J drw{r — fj) < ( — - — j > w{f — fj) is the onsite kinetic energy of the atoms. Pq = J dfw{f — 

fj) cos^ {k' x)w{f — fj) is the effective onsite potential energy of the atoms. The last term in the above Hamiltonian 

represents the two-body atom-atom coupling where v = ^^^^^ J dr\w{'f)\^ represents the effective onsite atom-atom 
interaction energy. 

III. LINEARIZATION OF QUANTUM LANGEVIN EQUATIONS 

Now, we study the quantum langcvin equations(QLEs) of the system. The QLEs of motion for the boson field 
operator aj, cavity photons b and movable mirror (mode) operator c are given as: 

ajii) = -i^aj{t) - i^{hvob^t)b{t) + Vci)aj{t) - i'^a\{t)aj{t)aj{t), (4) 
an n ■' 

b{t) = -iPoVob{t) a]it)aj{t) + iAcb{t) + ieiJmb{t)[c{t) + ct(f)] jj' - '^b{t) + y/Kbin{t), (5) 
j 

c{t) = -iwmc{t) + iewmb\t)b{t) - Tmc{t) + ^/vZU{t)■ (6) 

Due to leakage of photons from the mirror, the cavity field is damped. Therefore, we have introduced k as the 
cavity field damping rate. Also, the interaction of the cantilever with environment damps the mechanical mode with 
damping rate Vm- During the experiment, the condensate atoms are robust and therefore the depletion of atoms 
is not significant. Since the condensate temperature is much smaller than hn/kB, so, there is negligible coherent 
amplification of the condensate motion. The mechanical mode is also affected by a random Brownian force which has 
^ as zero mean value. The vaccum radiation input noise is represented by bin{t). The correlation functions for the 
input noise operators are given in the Appendix 1. The QLEs are linearized around the steady state as b{t) — > (3 + Sb{t), 

c{t) ^ 7+(5c(t) and a,(t) ^ ^^^'^ where (3 (= _,^^+^+,pJ^N-,2Re(^) '7 (^6(7) = ^|^) and are the 

steady state values of cavity mode, mechanical mode and condensate density respectively. N{= a}-{t)aj{t)) atoms 
occupy M number of lattice sites. It is assumed that all the sites of optical lattice are same, therefore aj{t) is replaced 
by a{t). By introducing the amplitude and phase quadratures as 5qa{t) = [5a{t) + (5a^(t)], Spa{t) = i[Sa^{t) — Sa{t)], 
5qt,(t) = [6b{t)+5b\t% 6pbit) = i[6b^{t)-db{t)], 5q{t) = [(5c(i)+fct(t)], 5p{t) = i[5c^ {t)-5c{t)], q^t) = [bin{t)+bl{t)] 
and Pin{t) = i[bl^{t) — bin{t)], the linearized QLEs for the system are given as follows: 

Sqait) = PiSpait), (7) 

dpa{t) = -02Sqait) - 2gMb{t), (8) 

^<ib{t) = -'^hb{t) + y/Rqin{t) - AdSpb{t), (9) 

SPbit) = -'^Spbit) - 2gMa{t) + 2Gpdq{t) + V^pin{t) + Addqt,{t), (10) 

5q{t)=^mSp{t), (11) 

6p{t) = -i^m6q{t) + 2Gp6qb{t) - r„5p(t) + W{t), (12) 
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where = -PovoN + A, + 2Gj, = PoVoPVN, G = ew„, ly = Ko/h + Povof3^ + PoV.i/H, U^ff = ^, 
Pi ~ v + Ueff, 1^2 ~ V + 'iUeff and W{t) = m,[^m{'^) ~ ^in(^)] which satisfies the correlation given in the Appendix 
I. 

Let us now understand how the radiation pressure cools the mechanical oscillator. Radiation pressure exerted 
on the movable mirror by the optical cavity mode forms a system which acts as another reservoir connected to the 
mechanical oscillator when the cavity is properly detuned. As a result, the effective temperature of the vibrational 
mode is the temperature between the two reservoirs i.e., between the initial bath temperature and temperature of this 
effective optical reservoir. This effective temperature is zero in practice. Hence the quantum ground state is acheived 
when the coupling to the initial reservoir, given by the damping rate of the movable mirror F™, is much smaller than 
the coupling to the effective optical reservoir. This explains us that the radiation pressure coupling should be strong 
for significant cooling of the mechanical oscillator. Now, we shall see in the next sections how the different cooling 
techniques help in cooling the mechanical oscillator to its quantum ground state. 



IV. BACK-ACTION COOLING 



Back-action dynamics have been realized in a diverse variety of physical systems, including those of ultracold atoms 
[31I [l^l-The randomness present in the unavoidable stochastic back-action forces is due to the photon shot noise. 
These forces arise from the radiation pressure. 

OPTICAL CAVrVY 

MOVING 

FIXED MIRROR 
MIRROR 



Figure 1: (color online) Setup for Back- Action Cooling. It configures an optomechanical system involving Bose Einstein 
Condensate (BEC) inside a Fabry-Perot cavity with one fixed mirror and another movable light-end mirror. The laser pump is 
used to drive the cavity mode. 

In this section, we evaluate effective frequency, effective damping rate and displacement spectrum for the mechanical 
oscillator (mirror) in the back-action cooling scheme. Also we show how the ground state cooling of the mirror is 
approached in this scheme. 

The displacement spectrum in fourier space is evaluated from 




5,(w) dLj'e-'^'^+'''^'{Sq{Lu)Sqiuj') + Sq{Lu')SqiLj)), (13) 

using the correlations in fourier space given in the Appendix I. Therefore, the displacement spectrum in the fourier 
space for the movable mirror is given as : 



Sgiuj) = \Xeff{uj)\^[SthH + Srp{uJ,Aa)], (14) 

where Sthi^) is the thermal noise spectrum arising from the Brownian motion of the mirror and Srp{uj, Ad) is 
the radiation pressure spectrum including the quantum fluctuations of the condensate. Xeffi^) is the effective 
susceptibility of the oscillator. 

Here 



Sth{u}) ^ — wcoth ( 77^^ ) , 
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where 



UJ h + (A^ + — - W j 



(16) 



(17) 



where 



(18) 



(19) 



Xe//(w) is the effective susceptibihty of the resonator ahered by the radiation pressure and condensate fluctuations 
with 



IXe//(a.)|' 



(20) 



The effective mechanical susceptibihty of the osciUator gives us the effective resonance frequency and effective 
damping rate as : 



where 



(21) 



_ AG^fj^AdO^rnji^' - PMico^ - /3i/32)(Ag + ^~oj^)^ 4ggAd/3i] 
""^ X{uj) 



(22) 



AG^(3^AdiOr,M^^ - 



(23) 



Eqn.(|2ipand (j22p show that the mechanical frequency of the mirror gets modified by the quantum fluctuations of 
mirror due to the radiation pressure and condensate fluctuations. This is the so-called optical spring effect. It can 
be observed that the frequency due to the optical spring term does not get altered significantly for high resonance 
frequencies, such as those of Refs.(l8l-l20j where ^ 1 MHz. 

Experimentalh^ the mirror may have mechanical frequency varying from 27r x 100 Hz [49|, 2tt x 10 kHz [35| o r 
2tt X 73.5 MHz [sOl with corresponding damping rate from 2tt x 10~^ Hz [11], 27r x 3.22 Hz [sl] or 27r x 1.3 kHz]50|. 
A high finesse Fabry Perot optical cavity having decay rate k — 271 x 8.75 kHz [5l| (27r x 0.66 MHz [s^]) consists of 
cloud of BEC with an order of 10^ ^"^Rb atoms [46| interacting with cavity field may have coherent coupling strength 
as go = 2n X 5.86 kHz [5l| {2n x 14.4 MHz [S^]). The loss of photons through the cavity mirrors decreases the energy 
of the cavity mode which further minimises the interaction of atom light field. This loss of photons can be reduced in 
high finesse optical cavities. 

In Fig. 2(a), the normalized effective mechanical frequency {uf^l^ /um) of the oscillating mirror in the absence of 
BEC (dashed line) is compared with the case in the presence of BEC (solid line) as a function of normalized frequency 
(oj/wm). This figure shows an extra resonance dip in the presence of BEC. Significant deviation from the bare 
frequency Wm is observed around uj — ±uJm- This deviation is enhanced in the presence of BEC. Fig. 2(b) illustrates 
the variation of normalized effective mechanical frequency of the mirror as a fuction of Lu/ujm for different two body 
atom-atom interactions, Ueff = 0.2a;m(dashed line), f/e// = 0.4ajm(dot dashed line) and Ue/f = 0.8aJm(solid line). 
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Figure 2: (color online) Plot of normalized effective mechanical frequency {uj^-^ /ujm) as a function of dimensionless frequency in 
back-action cooling scheme. Fig. 2(a) shows the variation of the dimensionless effective mechanical frequency in the absence of 
BEC (dashed line) and in the presence of BEC (solid line) with Ueff ~ O.ScJm- Fig. 2(b) represents the dimensionless effective 
mechanical frequency for three different values of atomic two-body interaction with Ueff ~ 0.2cjm (dashed line), Ueff = OAujrn 
(dot dashed line) and Ueff ~ O.Sojm (solid line). Other parameters used are : Tm = fO~^cJm, Ad — —ujm, n = O.ScJm, 
G = 4cj„,/3 = 0.05, ly = O.Oloj^.gc = O.SuJra and ksT/huj^ = 10^ 



a b 




Figure 3: (color online) Plot of normalized effective damping rate (rf,{-'/tJm) with dimensionless frequency in back-action 
cooling scheme. Fig. 3(a) represents the variation of dimensionless effective damping rate in the absence of BEC (dashed line) 
and in the presence of BEC (solid line) with Ueff = O.SoJm. Fig. 3(b) gives the deviation of dimensionless effective damping 
rate for three values of atomic two-body interaction with Ueff ~ Q.2u)m (dashed line), Ueff = QAojm (dot dashed line) and 
Ueff = O.SoJm (solid line). Other parameters chosen are the same as in figure 2. 

It can be seen from the figure that the deviation of the mechanical frequency of mirror from its resonance frequency 
Wm decreases with the increase in condensate two-body interaction since the condensate becomes more robust with a 
higher two-body interaction. With in a set of experimentally achievable parameters, the atomic two-body interaction 
can be modified using the condensate cloud having dimensions 3.3 /xm [Si] (290 nm [sHi) with length 20 ^m [sTJ (615.5 
nm [1^). It may also vary using scattering length ranging from lOao to 190ao (ao — Bohr radius) [52]. 

Fig. 3(a) shows the normalized effective mechanical damping [V^l^ /u)„i) of the mirror as a function of normalized 
frequency {uj/ujm) in the absence of BEC (dashed line) and the presence of BEC (solid line). As seen from this figure, 
below resonance (w < LOm), the effective damping of the oscillator is more in the presence of BEC whereas above 
resonance (w > Wm), it is less in the presence of BEC. Also, the plot of normalized effective mechanical damping 
of the mirror with normalized frequency is illustrated in the Fig. 3(b) for the three different values of condensate 
two-body interactions, Ueff = 0.2wm(dashed line), Ueff = 0.4wm(dot dashed Hne) and Ueff = 0.8a;m(sohd line). 
A higher condensate two-body interaction enhances the effective damping of the mirror below resonance whereas it 
decreases the effective damping of the mirror above resonance. It shows an exception for the case Ueff = QAuim 
as explained later. Fig. 4(a) represents the displacement spectrum Sq{uj) as a function of dimensionless frequency 
{uj/oJm) in the presence of BEC (solid line) and the absence of BEC (dashed line). We observe the normal mode 
splitting into two modes in the absence of BEC and we find that the normal mode splits up into three modes in the 
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Figure 4: (color online) Figure shows the displacement spectrum as a function of dimensionless frequency in back-action cooling 
scheme. Fig. 4(a) shows the displacement spectrum in the absence of BEC (dashed line) and in the presence of BEC (solid line) 
with Ueff = 0.8a;m. Fig. 4(b) represents the displacement spectrum for three different values of atomic two-body interaction 
with Ueff ~ 0.2ajm (dashed line), Ueff = OAujm (dot dashed line) and Ueff = O.Sujm (solid line). kBT/ftWm = 10'^ and the 
various other parameters used are same as in figure 2. 



presence of BEC. This extra mode is the result of additional quantum fluctuations of the condensate (Bogoliubov 
mode). Also below resonance, we observe that the amplitude of the peak(l) of displacement spectrum in the absence 
of BEC is greater than that in the presence of BEC. This is due to the higher damping rate in the presence of BEC as 
illustrated in Fig. 3(a). This conclusion is exactly opposite to that for the case above resonance. Fig. 4(a) depicts that 
the amplitude of peak (1) in the absence of BEC is more than with BEC whereas the amplitude of peak (2) in the 
absence of BEC is less than with BEC. This represents the energy exchange between different modes of the system. 
Also, Fig. 4(b) illustrates the displacement spectrum varying with normalized frequency for three different values of 
condensate two body interaction, Ueff = 0.2a;m(dashed line), Ueff = 0.4ajm(dot dashed line) and Ueff = 0.8a;m(solid 
line). Here coupling between the mechanical mode of the mirror, fluctuation of cavity field around steady state and the 
fiuctuations of the condensate (Bogoliubov mode) around the mean field results in normal-mode splitting (NMS). For 
the observation of NMS, it is significant to note that the time scale for the exchange of energy between the mechanical 
mode, photon mode and the Bogoliubov mode should be faster than the decoherence of each mode. The beat of 
laser pump photons with the photons scattered from the condensate atoms is the source of cavity field fiuctuations. 
Since the frequency of the Bogoliubov mode of the condensate is directly proportional to \/U7ff, fig-4(b) shows a 
variation in the displacement spectrum for different Uef / . The spectrum shifts towards the right with the increase in 
condensate two body interaction. Recently, an experiment reveals that Bogoliubov mode of cloud of ultracold atoms 
interacting with optical resonator has momentum ±2fcc(fcc is the cavity wave number) (3lj |. 

The condition r„i << uim << ksT/h is always considered in standard optomechanical experiments (2^ [53l - [55| . 
fn this limiting case, we consider the approximation that cot]i{huj/2kBT) ~ 2kBT/huj. Now, in order to acheive the 
ground state cooling of the mechanical resonator, we measure the average energy of the resonator in steady state, 
which is given by [56| 



[{Sq') + {5p')] ^ t^„, i^ueff + . (24) 

The system has to be stable in order to be in steady state. Hence, the stability conditions for back-action cooling 
scheme given in the Appendix H must always be satisfied in this regime such that all the poles of Xeffi^) lie in the 
lower complex half-plane. Here, equation (|24l) implies that the effective phonon number of the mirror's motion can 
be written as 



eff = 2 [i^l') + i^P') ~ 1] ' 



(25) 



where {5g^ and (Sp^) represent the displacement and momentum variances of the oscillator respectively which are 
defined as 
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Figure 5: (color online) Back-action cooling scheme: (a) Plot of the mirror's position variances (5?^) (Solid line) and the 
mirror's momentum variances (5j3^) (Dot dashed line) versus the normalized effective detuning (Ad/uim) in the absence of 
BEC and for three different values of atomic two-body interaction with Ueff ~ 0.4tJm, Ueff = O-SuJm and Ueff = t^m- Other 
parameters chosen are: Tm = W~^UJm, k = O.lwm, G = ujrn, P = 0.05, u — O.Olo;™, Qc ~ O.^LUm and kBT/HcOrn = 10^. (b) 
Effective phonon number with dimensionless effective detuning in the absence of BEC (thick Solid line) and for three different 
values of atomic two-body interaction with Ueff = OAujm (dashed line), Ueff ~ O.SiUm (thin Solid line) and Ueff = uum (dot 
dashed line). Other parameters chosen are the same as in (a). 




Figure 6: (color online) The effective phonon number versus the normalized atomic two-body interaction in back-action cooling 
scheme at Ad = —uum- Rest of the parameters are same as in figure 5. 



iSq')^ I (26) 

iSp') = r P^S,{.). (27) 

We have solved these oscillator variances numerically using the MATHEMATICA 8.0. Generally, there is no 
equipartition of energy, {Sq'^) 7^ (Sp^)- In order to approach the ground state cooling, the effective phonon number 
should be less than one i.e., n^ff < 1 which can only be attained if the initial average-thermal excitation number 

n = — 1 is not excessively large. This can be possible even at cryogenic temperatures if w,„ is 

sufficiently large. 

In order to get an intuitive picture, we have plotted both the displacement and momentum variances of the mirror 
as a function of dimensionless detuning Ad/w,„ in the absence and presence of BEC as illustrated in Fig. 5(a). Fig. 5(a) 
shows that, for large values of detuning, the variances have lesser values in the absence of BEC than in the presence of 
BEC. It is evident from the figure that one can access the gound state of the system by increasing the Bose- Einstein 
condensate two-body interaction. This implies that Ueff may alter the cooling process significantly. So, in order 
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Figure 7: (color online) The normalized effective mechanical damping rate versus the normalized atomic two-body interaction 
at the resonance frequency {uj = iOm) in back-action cooling scheme. Other parameters chosen are same as figure 2. 

to acheive the best condition for the ground state cooling of the system, Uef / can be used as an additional control 
parameter. Fig. 5(b) depicts the effective phonon number n^ff as a function of dimensionless detuning A^/wm for the 
absence of BEC (thick solid line) and the three different values of Ueff, Ueff — 0.4a;m(dashed line), Ueff — 0.8a;m(thin 
solid line) and Ueff = lwm(dot dashed line). Clearly, with the increase in Ueff, the minimum value oiueff is shifting 
towards zero such that, for Ueff = iojm, we have rieff — 0.52 at A^; = — l.lSwm- However, the least value of Ueff 
is reached in the absence of BEC which is nearly 0.51 at A^ — —lujm- Hence, in back-action cooling scheme, we 
have noticed that the minimum value of effective phonon number is found in the absence of BEC. Moreover, we have 
observed that for small values of A^;, presence of BEC gives better results for higher Ueff[= O.So;™ and ujm) than the 
absence of BEC. One can also observe from Fig. 5(b), the effective temperature of the oscillator does not vary much 
in the presence of BEC(C/e// — O.So;™ and oJm) as compared to that in the absence of BEC. This reveals that for 
a wide range of detuning(Ac;), the effective temperature of the mirror is not changing significantly for higher Ueff- 
Therefore, one can use condensate two-body interaction as a new handle to acheive and sustain a low temperature of 
the mirror with BEC over a wide range of A^. The atomic two-body interaction is proportional to number of atoms 
{N). In the recent past, [46| has shown experimentally that increasing the number of atoms enhances the damping of 
the oscillating membrane which is coupled to BEC through cavity field. This validates our work that the energy of 
the mirror decreases by increasing Ueff- Further, it is noticed that both the variances tend to {5q^) — {Sp^) ~ 1/2 
for G — 0.5wm (keeping other parameters same) i.e., energy equipartiton is satisfied in this chosen parameter regime 
but it results in the approximate same ground state oscillator energy as we have acheived in the above mentioned 
general case of no energy equipartition {(5q^^ ^ (^Sp^'j)- We find that the best cooling regime is acheived in the good 
cavity limit condition (k << uJm) by taking k — O.lwm as mentioned in 12]. 

Variation of effective phonon number with Ueff/i^m at A^ = —uJm is also shown in Fig. 6. It depicts a rapid 
increase in Uef / for a range of Uef / . This corresponds to the sudden decrease in the effective damping rate of mirror 
as illustrated in Fig. 7. It is by virtue of the fact that, in the range < Ueff < 0.57a;m, the second term in Eqn. 
increases as (w^ — /3i/32) decreases. On the other hand, for Ueff > 0-57u!m, this term decreases since (o;^ — /3i/32) 
becomes negative. It explains the exception observed in Fig. 3(b) and 4(b) for Ueff = 0.4wm. In the next section, 
we study the system using the cold damping feedback technique. 



Cold damping feedback technique is an alternative method to improve the cooling of a mechanical oscillator b 



displacement of the oscillator is measured through phase-sensitive homodyne detection of the cavity output which is 
fed back to the resonator using a force proportional to the oscillator velocity [13, 1131 ■ 
Therefore, the QLEs involving a feedback force for this scheme are given by : 



V. 



COLD DAMPING FEEDBACK SCHEME 




6qa{t) = l3lSpa{t), 



(28) 



6pa{t) ^ ~P2Sqa{t) ~ 2ge5qb{t), 



(29) 
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FEEDBACK 

Figure 8: (color online) Setup for Cold Damping. The configuration is same as in case of Back-Action Cooling involving an 
additional feedback loop with a force proportional to oscillator velocity. Cavity output field is homodyne detected using the 
beam splitter. 



Sqb{t) = --Sqb{t) + y/Kqin{t) - AdSpb{t), 



(30) 



Spb{t) = -^^Vb{t) - 2gMa{t) + 2G(36q{t) + VKPm(t) + AdSqbit), 



Sq{t) = uj„iSp{t), 



(31) 
(32) 



Sp{t) - -UJmSq{t) + 2Gl3Sqb{t) - TrnSp{t) + W (t) 



dsg{t - s)Spest{s)- 



(33) 



where Pi = + Ueff , 1^2 ^ v + SUeff and the filter function g{t) is the causal kernel chosen as [12| 



(34) 



such that 



5(w) 



I — iuj/ujfb^ 



(35) 



where g{uj) is the Fourier transform of g{t). Here gcd is the positive feedback gain and A quantifies the homodyne 
detection efficiency of the photodetector. Generally the detector efficiency is A < 1 for the inefficient homodyne 
detection if the effect of additional noise is considered, but for ideal homodyne detection, A = 1 [s^l • ojji^ signifies the 
feedback loop delay time. Moreover, Spest{s) represents the estimated intracavity phase quadrature which is evaluated 
by using the generalized input-output phase relation given by [57l - l59j 



ut(i) = \/n\/\5pb{t) - VXpin{t) - \/l - Xpvit). 



(36) 
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Here Pout (t) and pin (t) represent the output homodyne phase quadrature and input noise respectively. Also py (t) 
is the associated vaccum field quadrature which describes the additional noise for A < 1 such that 



Spest{t) 



Poutjt) _ ^gp^^^-^ _ ^Pinjt) _ Vl - \pv{t) 



(37) 



where Pv{t) and Pin(t) are uncorrelated. Now the QLEs supplemented with the feedback term are solved in the 
frequency domain such that the displacement and momentum variances of the oscillator are given by the Eqns.p6|) 
and (|27p respectively. It gives a distinct displacement spectrum using the correlations given in Appendix I. Explicitly, 
it can be written as 



(38) 



Here, Sth{^) and Srp{uj, A^) are given by the Eqns. (llSp and (jl6p respectively. In this scheme, the position spectrum 
consists of an additional feedback-induced term expressed as 



Sfb{uj) 



lGp^{^^ - p,h) {a,(c.2 - ^i/?2)(w2 _ - A2) + 4g2/3i(A2 + 4)} {g{u:) + g{-o,)} 



(39) 



Xiu;) 



X{u) 



which arises since the cold damping loop feeds back the measurement noise into the dynamics of the movable mirror. 
Xg^j(aj) is the effective mechanical susceptibility of the movable mirror, modified by the filter function, given by 



xf//(-) 



LOr, 



where 



[«-c^2 + i^r„,) + xf(^)]' 

2G/?w„(w2 _ p^p^-^ |g(-c^)^/A(^t^ + f ) + 2G/3Ad} 
It gives us the effective resonance frequency and damping rate using Eqn. ([35]): 



(40) 



(41) 



^ni (W) = [W„j + W„f J 

where 



op,cd 



Xiiuj) 



+ Xi(w)(w2 



■ p\P2) I n , — T^y^ 



2uj'^gcdi^fbX K 



(^2+^2 ) ^2 



(42) 



(^2-/3i/32)[A2 + --c.2]-4g,^A,/3i 



'^gcd^fbX ^ 2 I ^/fc^^ 

2 



1/32) ( 



Xi(cj)At(w^-/3i/32) \^GIiAd + 



2Lj'^gcdt^fbX ,K 



fb) : 



(43) 
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Figure 9: (color online) Plot of normalized effective mechanical frequency [oj'^^ /uj^n) as a function of dimensionless frequency in 
cold damping feedback scheme. Fig. 9(a) shows the variation of the dimensionless effective mechanical frequency in the absence 
of BEC (dashed line) and in the presence of BEC (solid line) with Ueff ~ O.Sujm- Fig. 9(b) represents the dimensionless effective 
mechanical frequency for three different values of atomic two-body interaction with Usff = 0.2u}m (dashed line), Ueff ~ OAuurn 
(dot dashed line) and Ueff = O.Sujm (solid line). Other parameters used are : Tm = 10~^uim, Ad = —uJm, n — 1.2wm, G = GoJm, 
P = 0.05, V = 0.015a;™, = 0.3tJ,„, cj/t = 4cj™, = 0.8, A = 0.8 and ksT/hwrn = 10^ 



a b 




0.8 0.9 1.0 1.1 1.2 1.3 0.8 0.9 1.0 1.1 1.2 1.3 



Figure 10: (color online) Plot of normalized effective damping rate {V^^ /ujrn) with dimensionless frequency in cold damping 
feedback scheme. Fig. 10(a) represents the variation of dimensionless effective mechanical frequency in the absence of BEC 
(dashed line) and in the presence of BEC (solid line) with U^ff ~ O.Sujm- Fig 10(b) gives the deviation of dimensionless effective 
mechanical frequency for three values of atomic two-body interaction with Ueff ~ 0.20;™ (dashed line), Ueff = OAujm (dot 
dashed line) and Ueff ~ O.ScJm (solid line). Other parameters chosen are same as in figure 9. 



where 

X{uj) 

In order to overdamp the mechanical mirror oscillations, an additional viscous force known as the feedback force is 
applied in this technique which is possible only when the estimated intracavity phase quadrature Spe^t is proportional 
to the oscillator position Sq{t). It can be achieved in the bad cavity limit where k > ujm- Hence, for cold damping 
feedback, we limit our discussion to the bad cavity limit. We have plotted the effective resonance frequency [uj'^^''^'^) 
and damping rate (r^-^''^^) as a function of dimensionless frequency (uj/uJm) in Figs. (9) and (10) respectively in order 
to compare them with the back-action cooling plots. It can be seen again from the fig. (9) that there is no significant 
shift in the frequency in both absence and presence of BEC in the chosen parameter regime. Hence, there is negligible 
optical spring effect as dominates over uj"^''^'^ for higher resonance frequencies. However, the mechanical damping 
rate shows significant variation with change in frequency. Below resonance, effective damping increases by adding 
BEC to the system [see Fig. 10(a)] and by increasing the atom-atom interaction [see Fig. 10(b)] with an exception for 
Ueff = OAujm- This exception in the result can be explained from the Fig. (11) which clearly shows the variation 
in mechanical damping rate with increasing Ueff. The reason for this sudden decrease in the damping rate in a 
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Figure 11: (color online) The normalized effective mechanical damping rate versus the normalized atomic two-body interaction 
at the resonance frequency {uj — ujm) in cold damping feedback scheme. Other parameters chosen are same as figure 9. 
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Figure 12: (color online) Figure shows the displacement spectrum as a function of dimensionless frequency in cold damping 
feedback scheme. Fig. 11(a) shows the displacement spectrum in the absence of BEG (dashed line) and in the presence of BEG 
(solid line) with f/e// = O.Suum- Fig. 11(b) represents the displacement spectrum for three different values of atomic two-body 
interaction with Ueff = 0.2ujm (dashed line), Ueff = OAum (dot dashed line) and Ueff = 0.8u!m (solid line). kBT/hujm = 10"^ 
and the various other parameters used are same as figure 9. 



particular region of Ueff for cj — ujm is similar as we have described in case of back-action cooling. However, in the 
zero detuning case, one can clearly observe from the expressions ([42]) and ([43]) that uj'^l^''^'^ and Tff^-'^'^ are independent 
of BEC parameters and behave in the same manner as in the case for absence of BEC. In the adiabatic limit for zero 
detuning i.e., k, w/f, >> w, we get u)^^''^'^ ~ and V^^''^'^ = r„ -I- 9cdG^^hrA^ This shows that the effective damping 
rate increases in this scheme without involving any significant change in the resonance frequency of the oscillator. 
We have also shown the plot of displacement spectrum {S'^'^{uj)) as a function of dimensionless frequency (oj/uJm) in 
Fig. 12 to compare it with the corresponding curve for the back-action cooling. This scheme involves an additional 
feedback-induced term denoted by ^/^(a;) which gives a distinct displacement spectrum than the back-action cooling 
scheme. This feedback induced term Sfb{i-^) can be used as an additional handle in the cold damping feedback scheme. 
The displacement spectrum in this scheme can be manipulated by a coherent control over the feedback parameters gcd 
and ujfb- The variation in amplitude of the peaks in the absence and presence of BEC of the displacement spectrum, 
shown in fig. 12(a), represents the energy exchange between the different modes of the system. Only two-mode splitting 
is observed in the absence of BEC [see Fig. 12(a)] while normal mode splits up into three modes due to atom-atom 
interaction [see Fig. 12(b)]. 

Now, we characterize the steady state energy of the mirror in the cold damping feedback scheme. The optimal 
cooling conditions can be obtained using the stability conditions which are modified in this scheme given in the 
Appendix II. The oscillator variances represented by Eqns. (|26p and ([27]) and the mean energy of the oscillator given 
by Eqn. (|M]) . using the corresponding position spectrum for the cold damping case, are calculated numerically with 
the help of MATHEMATICA 8.0. This shows that, also using cold damping, there is no equipartition of energy 
i.e., (Sq^) (Sp^) in general. Fig. 13(a) illustrates both the variances {Sq^) and {Sp^) as a function of dimensionless 
detuning (A^/wm) for the absence of BEC and three different values of condensate two-body interaction (Ueff), 
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Figure 13: (color online) Cold damping feedback scheme: (a) Plot of the mirror's position variances {Sq^'j (solid line) and 
the mirror's momentum variances {Sp^) (dot dashed line) versus the normalized effective detuning (Ad/uim) in the absence of 
BEC and for three different values of atomic two-body interaction with Ueff ~ OAidm, Ueff = O-Sujm and Ueff = ijJm- Other 
parameters chosen are: Fm — 10~^a;m, n = 5^;^, G = 2uJm, /3 = 0.05, u = ujrn, Qc = oJm, cofb = 5u!rn, Qcd = 0.4, A — 0.85 and 
ksT /hujm = 10^. (b) Plot of effective phonon number with dimensionless effective detuning in the absence of BEC (thick Solid 
line) and for three different values of atomic two-body interaction with Ueff ~ OAujm (dashed line), Ueff ~ O.Sujm (thin Solid 
line) and Ueff = (jJm (dot dashed line). Other parameters used are same as in fig(a). 
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Figure 14: (color online) The effective phonon number versus the normalized atomic two-body interaction at — —ujm in 
cold damping feedback scheme. Rest of the parameters are same as figure 13(a). 



Ueff = OAuJm, Ueff = O.Swm End Ueff = Iwm- As It Can be seen from the figure, the displacement and momentum 
variances decrease significantly on adding BEC to the system. Moreover, both variances decrease with the increase 
in condensate two-body interaction. However, as approaches zero, oscillator variances increase drastically in the 
presence of BEC. Therefore, one can infer that the presence of BEC helps in cooling the mirror for smaller detunings. 
So, in order to approach the optimal cooling conditions in cold damping scheme, the sum of these variances given by 
equations (|26|) and ()27|) is to be minimized which can be done by increasing the condensate two-body interaction. The 
effective phonon number given by Eqn. (j25p is shown in Fig. 13(b) as a function of dimensionless detuning (A^/wm) for 
the absence of BEC (thick solid line) and the three different values of atom-atom interaction, Ueff = 0.4wm (dashed 
line), Ueff = 0.8aJm(thin solid line) and Ueff = lWm(dot dashed line). From the figure, we observe that the minimum 
value of effective mean excitation number is reached for Ueff = 0.8cj„ which is 0.658 with A^; = — 1.5wm. Although, 
as Ac; -> 0, the cavity field amplitude quadrature dqb{ui) becomes insensitive to the mirror motion and ^e// increases 
considerably with increasing Ueff- So, the ground state cooling can be approached by adding BEC to the system 
involving cold damping feedback scheme for lower detunings. Hence, by varying the atom-atom interaction, one can 
optimize the cooling process to acheive the ground state. Using cavity self-cooling scheme, least value of Ueff is 
acheived in the absence of BEC while using cold damping feedback technique, minimum value is obtained by adding 
BEC to the system. Moreover, with BEC, we are getting the least value of effective phonon number using back-action 
cooling scheme. We have also shown the variation of effective phonon number rieff as a function of Ueff/uJm in 
Fig. (14). It clearly illustrates the fact that, in the presence of BEC, the ground state cooling of the mechanical 
oscillator can be achieved for higher condensate two-body interaction only. This sudden increase in the value of Uef / 
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for very small atom-atom interaction is the consequence of the decrease in the effective mechanical damping rate of the 
mirror given by Eqn. (j43p in this region. Since in this chosen parameter regime, the second term in Ean. (l43l) becomes 
excessively small and even negative for very small Ueff- Moreover, we have also examined that by taking ujfb — 0.5aJm, 
keeping other parameters same as before, the uncertainty principle is stratified for both displacement and momentum 
variances. This implies that {5q^) ~ {Sp^) — 1/2, therefore, energy equipartiton is satisfied. However, the ground 
state could not be approached in this regime i.e., we are getting rieff > 1. 

VI. CONCLUSION 

In this paper, we have studied how the back-action cooling and cold damping feedback schemes help in cooling the 
mirror to its quantum ground state by using a BEC confined in an optical cavity. The atom-atom two body interaction 
can be used as a new handle to cool the quantum device in both the schemes. It provides a systematic control of the 
system which can be altered either by using number of atoms or s-wave scattering length. Both the techniques show 
distinct displacement spectrums involving energy exchange between the different modes of the system. A coherent 
control over the feedback parameters, Qcd and cj/f,, can manipulate the displacement spectrum in the cold damping 
feedback scheme. In back-action cooling, the least value of effective phonon number is obtained without BEC whereas 
in cold damping, it is obtained with BEC. The effective temperature of the oscillator does not vary much for higher 
condensate two-body interaction with self-cooling. Also in this scheme, higher negative detuning gives lesser value of 
effective phonon number in the presence of BEC as compared to that in the absence of BEC. As we compare both 
the schemes, the minimum value of effective phonon number is obtained using cavity self-cooling with BEC. We have 
analyzed that both the techniques help in approaching the quantum ground state of the oscillator. Dicrimination in 
the ideal cooling conditions for these schemes exhibits that back-action dynamics is more convenient in good cavity 
limit (k << w„j) while cold damping is preferable in bad cavity limit (k >> uj.„i). 
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VIII. APPENDIX I 



Following correlations are satisfied by the input noise operators [10|, [161, [571, [6 



{h^n(t)b,r^{t')) = {bl{t)hn{t')) = 0, 



(45) 



(46) 



The noise operator due to the Brownian motion of the mirror is given as W{t) = i\/^m[^m{^) ~^m{t)] which satisfies 
the following correlation [57| : 



{W{t)W{t')} = ^ [ |^e-'"(*-*') 



1 + coth 



\2kBT 



(47) 



where T is the finite temperature of the bath connected to the cantilever and ks is the Boltzmann constant. 
Brownian noise is the random thermal noise which arises from the stochastic motion of the mechanical oscillator 
(mirror). It should be noted that Brownian noise is non-Markovian in nature. The quantum Brownian motion of the 
movable mirror give s rise to the thermal noise term in the measured phase noise spectrum of the optical field reflected 
from the cavity [lOl IT6| . 

The amplitude and phase quadratures of the input noise operator satisfy the following correlations in Fourier space 

M- 

(gm(w)ftn(a;')) = 2TrS{uj + uj'), (pi„(tj)pi„(a;')) = 2nS{uj + uj'), {qi„{uj)pi„{uj')) = 2171 S {uj + uj'), {pin{uj)qin{uj')) = 
-2iTT5{uj + uj'). 
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Also in fourier space, the correlation function for the brownian noise operator is given as [l^l : 
{Wiuj)W{u;')) = 2tt^uj 1 + coth (5^) '^('^ + w') 

The correlation of vacuum field quadrature in the Fourier space is {p^,{Lo)py{Lj')) = 2t:5{ll! + lj'). 

IX. APPENDIX II 

In the back-action cooling scheme, two non-trivial stability conditions for the system in the absence of BEG, obtained 
by applying the Routh-Hurwitz criterion, are given as: 

Si=ao> 0, (48) 
52 = (030201 — 040^ — 0003) > 0, (49) 



where 



ao =^™(^ + A2)+c^„/32G2Ad, (50) 



^.2 

01 = r™(A2 + — ) + w^K, (51) 



02 = — + A2+r,„« + c^^^, (52) 

03 = r,„ + K, (53) 

04 = 1. (54) 

While the stability conditions, evaluated by applying the Routh-Hurwitz criterion, for the system in the presence of 
BEG are as follows: 

53 = 60 > 0, (55) 

Sa = (65^463 + ^6^165 - b^hl - b^bl) > 0, (56) 



where 



60 = ^'^"f"^" + ^^/3i/32A^ + 4G2/32c,,„Arf/3i/32 + ^gl^aPiu^l, (57) 



&i = /3i/32KW™ + ^'^'f ^"^ + A/32r™A^ + Agl^dPiTn,, (58) 



h2 = ^ + c.2„A^ + 4G2/32c,„Ad + /3i/32r™« + + P.h^l + ^gl^dPi + fiiP2^l, (59) 



4 



^2 

63 = /3i/32(r™ + + r,„(A2 + _) + Kc^^, (60) 



^2 

64 = u;^-t-r™K+ — + A2+/3i/32, (61) 
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6r = r„, + K, 



(62) 



be = l. (63) 

Moreover, the Routh-Hurwitz criterion is equivalent to the condition by imposing all the poles of mechanical suscep- 
tibility ix'e'ffi^)) lower complex half plane. Hence, one gets a non-trivial modified stability condition for the cold 
damping case in the presence of BEC as follows: 



where 



S5 = (C6C3 + C2C5 - C5C4C3 - C6C5C1) > 0, 



(64) 



ci = -«(A^ + —){ujfbuj'^^ + a;/fc/3i/?2 + r„^i/32) - i/3i/32('^/6W^i + T„,KUJfb + ujI^k) ^^^^ 
-i-ig^Adl3i{ujfb + Tm) ~ i2Gl3u!rngcdl3il32Xujfb - 140^/3^ ArfW/bW„, 



C2 = (A^ + + TmCJ/b -I- I3il32) + /3l/32(w^ + FrnW/b -I- W/fcK -|- TmH) 



(66) 



C3 = z 



(A^ + — )(a;/fc + r,„) + LOfb{i^m + r^K + /3i/32) + Ka;,^„ + A^2(k + r„J + 2GPu>mgcd\i^fb 



(67) 



C4 



(A^ + y) + + t^/b(K + r„) + + uTr, 



(68) 



C5 = -i [w/b -I- K -I- Tm] , 



(69) 



Cg = 1. 



(70) 



In the absence of BEC, the stability condition for the system involving cold damping feedback can be obtained by 
putting all the BEC parameters to be zero in Eqns. (I64|) - (l70l) . 
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